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Abstract 

The phonon and thermodynamic properties of rare-earth-aluminum intermetallic s AIRE (RE=Y, Gd, Pr, 
Yb) with B2-type structure are investigated by performing density functional theory and density functional 
perturbation theory within the quasiharmonic approximation. The phonon spectra and phonon density of 
states, including the phonon partial density of states and total density of states, have been discussed. Our 
results demonstrate that the density of states is mostly composed of Al states at the high frequency. The 
temperature dependence of various quantities such as the thermal expansions, the heat capacities at constant 
volume and constant pressure, the isothermal bulk modulus, and the entropy are obtained. The electronic 
contribution to the specific heat is discussed, and the presented results show that the thermal electronic 
excitation affecting the thermal properties is inessential. 
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1 Introduction 

Aluminum alloys are widely used for aerospace industries, aircraft automotive industries, electronic industries and 
buildings due to light weight, good corrosion resistance, reasonably high strength and favorable economics. Addition 
of rare earth (RE) elements to Al-based metal alloys have received great attention since it can result in increasing their 
wear resistance |T], mechanical properties |2], electrochemical behavior (3) and thermal stability [4], etc. Thus, they are 

important for design of novel materials and further scientific and technical investigations. The development and design 
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of new Aluminum alloys need more fundamental physical data of these alloys. Recently, some investigations focusing on 
AIRE, AI2RE, and AI3RE systems intensely are performed from first-principle calculations based on density functional 
theory (DFT) E [3 El [9] [lOl [TH [T21|. The B2-type A1Y is a high-temperature stable while the temperature higher 
than 1473K (T3ll . Such metastale phases are also available for AlGd, AlPr, and AlYb, etc. intermetallics with B2-type 
structure lfT4l[T5l[T6l . Srivastava et. al. (5) investigated the electronic and thermal properties by using the tight binding 
linear muffin tin orbital (TB-LMTO). Tao et. al. (§1 [7) performed the projector augmented wave (PAW) method to 
calculate the elastic properties and predicted thermodynamical properties of B 2- AIRE intermetallics below 400^ within 
the Debye-type model. However, the accurate calculations of thermodynamic properties for B2-type AIRE intermetallics 
from the lattice dynamics, i.e. exact phonon spectra, have received little attention. 

Thermodynamical properties, such as entropy, specific heat, thermal expansion, and temperature dependent equa- 
tion of state (EOS) of AIRE intermetallics are very important to investigate their properties as well as applications. A 
simplified method for thermal expansion calculations is in the framework of density functional theory (DFT) with a 
Debye-Grimeisen based model fTTlL which is based on the long- wave approximation. This model is useful for calculating 
the thermal properties at low temperature, and thermodynamic properties of B2-type AIRE intermetallics have been car- 
ried out based on this model by Tao et. al |7 |. A more accurate approach has been made possible by the achievements of 
density functional perturbation theory (DFPT) from first principles |[T8lL which allowed exact calculations of vibrational 
frequencies in every point of the Brillouin Zone [19]. The vibrational free energy can be obtained using the quasiharmonic 
approximation (QHA). Moreover, QHA method lets one take into account the anharmonicity of the potential at the first 
order: vibrational properties can be understood in terms of the excitation of the noninteracting phonon. QHA based on 
DFPT provided a reasonable description of the thermodynamical properties of many bulk materials below the melting 
point JT9l [20j EU [22l [23]|, and had been applied with great success to more and more complex materials such as alloys 
[N1AI3 l24l l. perovskite [MgSiC>3 l25ll 1. and hexaborides [LaB6 and CeB6 (26)]. More recently, we have successfully 
investigated the thermodynamical properties of NiAl, YAg, and YCu l27l . and MgRE intermetallics |28 ] with B2-type 
structures. 

In this paper, we apply first principles calculations within QHA based on DFPT to study the thermodynamic properties 
of AIRE (RE=Y, Gd, Pr, Yb) intermetallics with B2-type structures. The phonon spectra and phonon density of states, 
including the phonon partial density of states and total density of states, have been discussed. Thermal expansions, 
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temperature dependence of isothermal bulk modules, heat capacities at constant volume and constant pressure, and the 
entropy as a function of temperature are presented. 

2 Theory 

To study the effects of changing temperature, one has to look at the Helmholtz free energy, incorporating the effects 
of thermal-electronic excitations and thermal- vibrations (phonons). The Helmholtz free energy at temperature T and 
constant volume V is give by (29l [30) 

F(V, T) = E (V) + F Q i(V, T) + F vib (V, T\ (1) 

where Eq(V) is the static contribution to the internal energy at volume V and can be easily obtained form standard DFT 
calculations. F e \(V, T) given in Eq. ([I]) is the thermal electronic contribution to free energy and is given by F e \ = E e \-TS e \. 
Here, The electronic excitation energy F el is given by 



r*oo r*£ F 

E Qi (V,T)= n(e,V)f(s)sds- n(s,V)sds, 
Jo Jo 



(2) 



where n(s, V) is the electronic density of state (DOS) at energy s and volume V, f(s) is the Fermi-Dirac distribution 
function, and sp represents the Fermi level. The electronic entropy 5 e i is formulated as 



S d (V, T) = -k B f n(s, V)[f(s)lnf(s) + (1 - /(e))ln(l - f(e))]de, (3) 
Jo 

where is the Boltzmann constant. Thermal electronic contribution to free energy is generally considered to be negligible 
away from the melting point of the material under consideration. 

F v fo(V, T) Eq. ([T]) is the vibrational free energy which comes from the phonon contribution. Within the quasiharmonic 
approximation (QHA), F^CK T) is given by 

F vib (K T) = k B T g In |2 sinh (^g^)j. (4) 

Here, the sum is over all phonon branches A and over all wave vectors q in the first Brillouin zone, % is the reduced Planck 
constant, and co q A(V) is the frequency of the phonon with wave vector q and polarization A, evaluated at constant volume 
V. 



The vibrational specific heat Cy at constant volume in the QHA from the following equation 

fcosh^^ 2 



c v =ZM^-) C0Sh (^-)- (5) 

qA 



The electronic specific heat can be obtained from 



^el 



and we denote that total specific heat at constant volume is then CV = C y h + Cy 1 . Due to anharmonicity, the specific 
heat at a constant pressure, C p , is different from the specific heat at a constant volume, C v goes to a constant which 
is given by classical equipartition law: C v = 3Nk B , where N is the number of atoms in the system while C p , which 
is what experiments determine directly, is proportional to T. QHA lets one take into account the anharmonicity of the 
potential at first order: vibrational properties can be understood in terms of the excitation of the noninteracting phonon. 
The equilibrium volume at temperature T is obtained by minimizing Helmholtz F with respect to V, i.e., miny[F(V r , T)]. 
The volume thermal expansion coefficient is given by 

a(T) = l v(%) P > (7) 

and the linear thermal expansion is described as 

a(T) - a(T c ) 

e(T) = — , (8) 

a(T c ) 

where a(T c ) is equilibrium lattice constant a(T) = [V(T)] 1/3 at T c = 300K. 
Then, C p can be obtained from Cy and a by 

C p = C v + a 2 BVT, (9) 

where B(T) = -l/Vd 2 F/dV 2 is the bulk modulus. 

The vibrational contribution to the entropy of the crystal is given by 

fico qA (V)\ hco qA (V) nco qA (V) 



Svib = - kB y m ( 2 sinh - ^2 coth 

^ \ 2k B T ) 2k B T 2k B T 

qA 



(10) 



3 Computational details 

In present work, the static energy and the thermal electronic contribution to the Helmhotz free energy are computed 
by using the first-principles calculations in the framework of the density-functional theory (DFT). We employed the gen- 
eralized gradient approximation (GGA) in the Perdew-Burke-Ernzerhof (PBE) 1 3 1 . 32] exchange-correlation functional 
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as implemented in the VASP code (33l|34j|35l. The ion-electron interaction is described by the full potential frozen-core 
projector augmented wave (PAW) method (36]|33, with energy cutoff of 600eV for plane waves. The Brillouin zones of 
the unit cells are represented by Monkhorst-Pack special k-point scheme (38) . We used a sample as 21 x21 x21 &-point 
mesh in the full Brillouin zone giving 726 irreducible ^-points to calculate the initial structures and electronic DOS. The 
radial cutoffs of the PAW potentials of Al, Y, Gd, Pr and Yb were 1.40, 1.81, 1.58, 1.64, and 1.74 A, respectively. The 3s 
and 3p electrons for Al, the 4s, 4p, Ad and 5 s electrons for Y, the 4f and 6s electrons for Gd, Pr, and Yb were treated as 
valence and the remaining electrons were kept frozen. The thermal electronic energies and entropies are evaluated using 
one-dimensional integrations from the self-consistent DFT calculations of electronic DOS using Fermi-Dirac smearing as 
shown in Eqs. ([2]) and ([3]). In order to deal with the possible convergence problems for metals, a smearing technique is 
employed using the Methfessel-Paxton scheme |39], with a smearing with of 0.05eV. 

The vibrational free energy were obtained from the first-principles phonon calculations by using PHONOPY (23j 
l40l [4TTl which can support VASP interface to calculate force constants directly in the framework of density-functional 
perturbation theory (DFPT) l42l . Phonon calculations were performed by the supercell approach. Since the chosen 
supercell size strongly influences on the thermal properties, we compare the vibrational free energies of 3 x 3 x 3 supercell 
with those of 5 x 5 x 5 supercell at 300K and 1000K, and find that the energy fluctuations between 3x3x3 and 5x5x5 
supercells are less than 0.01%. Hence, we chose the 3 x 3 x 3 supercell with 54 atoms to calculate phonon dispersions. 
We carried out DFPT calculations on this 54 atoms supercell using PBE-GGA exchange-correlations effects and 7x7x7 
&-point grid meshes for Brillouin zone integrations. 

In order to get the equilibrium lattice volume as a function of temperature, we have calculated total free energy at 
temperature points with a step of IK from to 1200K at 13 volume points. At each temperature point, the equilibrium 
volume V(T) and isothermal bulk moduli B(T) are obtained by minimizing free-energy with respect to V from fitting the 
integral form of the Vinet equation of state (EOS) | 43 ] at p = 0. These procedures applied for A1Y are demonstrated in 
Figure [T] where the Helmholtz free energy F(T, V) as a function of unit-cell volume at temperatures are shown. While, 
C p has be calculated by polynomial fittings for Cy and by numerical differentiation for dV/dT to obtain a(T). 
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4 Results and discussions 

4.1 Phonon spectra and density of states 

We calculated phonon frequency using the DFPT with force-constant method l44l with forces calculated using VASP. 
The calculated phonon spectra of A1Y, AlGd, AlPr, and AlYb, which are calculated by using 3x3x3 supercells, are 
displayed in Figures [2] Due to the crystal symmetry, the spectra curves are shown along high-symmetry direction G - 
X-M-T-R of the Brillouin zone. These dispersion curves have common framework. Since the B2-type intermetallics 
contain two atoms per primitive cubic unit cell, there are three acoustical branches and three optical branches. Our results 
obviously show that the degenerating of two transverse-acoustical (TA) modes and two transverse-optical (TO) modes in 
the G—X and G-R directions. At the G point, the degenerating of longitudinal-optical (LO) and TO modes is obtained, so 
there are no polarization effects in the B2-type AIRE intermetallics. The phonon DOS including the partial DOS (PDOS) 
and the total DOS (TDOS) are shown in Figures[3] Sampling a51x51x51 Monkhorst-Pack grid for phonon wave vectors 
q is found to be sufficient in order to get the mean relative error in each channel of phonon DOS. The phonon band gap of 
A1Y starts at 5 THz, while those of AlGd, AlPr, and AlYb starts at less than 4Hz. Among the four intermetallics, A1Y and 
AlYb has the greatest and lowest values of the maximum value of acoustic modes, respectively, since the atomic mass of 
Y is lightest and that of Yb is heaviest in the four calculated RE elements. The flat regions of phonon-dispersion curves, 
which correspond to the peaks in the phonon PDOS, indicate localization of the states, i.e., they behave like "atomic 
states" l23l . In the four intermetallics, the density of states are mostly composed of Al states above the phonon band gap 
since its atomic mass is much lighter than those of the rare earth elements RE (RE=Y, Dy, Pr, Tb). 

4.2 Bulk properties and thermal expansion 

The calculated results of the equilibrium lattice constants ao, the isothermal bulk modulus Bo, and the pressure 
derivatives of the isothermal bulk modulus B' at T = OK for AIRE (RE=Y, Dy, Pr, Tb) intermetallics together with 
the previous calculated results | 6] and the available experimental values (l3l[T5l[T6| are shown in Table [T] Our calculated 
results for the equilibrium lattice constants at T = OK show excellent agreements with the previous theoretical and 
experimental results. For the isothermal bulk modulus Bq and the pressure derivatives of the isothermal bulk modulus B' 
at T = 0, which are obtained from Vinet equation of state, the present results are within 1.7% errors from the previous 
calculated values obtained from the Rose's equation of state |6|. The temperature dependent isothermal bulk modulus 
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B(T) are shown in Figure]?] Among the four intermetallic compounds, through the temperature range - 1200^, lihgt 
RE A1Y and heavy RE AlYb have the highest and lowest bulk moduli, respectively, i.e, A1Y is the most incompressible 
and AlYb is the most compressible. With increasing temperature, the differences among the bulk moduli almost remain 
unchanged. The overall observation is that B T of the intermetallic s decrease with increasing temperature, and approach 
linearity at higher temperature and zero slope around zero temperature. 

As shown in Figure [T] the equilibrium volume at any temperature corresponds the minimum values of the fitted 
thermodynamic functions. The thermal expansion is observed as an increase in the equilibrium volume. The linear 
thermal expansion e defined by Eq. ([8} as a function of the four AIRE intermetallic s are shown in Figure [5ja). The linear 
expansions of the four compounds are found AlYb>AlY>AlGd>AlPr, and those of A1Y and AlGd are nearly equivalent. 
The coefficients of the volume thermal expansion a as a function of temperature are shown in Figure[5jb). With increasing 
temperature, the thermal expansion grows rapidly up to ~ 200^, and the slops become smaller and nearly constant at high 
temperatures except AlYb. In addition, a can be used to estimate the anharmonic effects according to Griineisen relation 
a = yC v /VB T , and Griineisen parameter y = -din oj/ din V. In QHA, the phonon frequencies at given lattice parameters 
are independent of temperature. In real crystal, it is not the case. The accuracy of QHA applied to AIRE (RE=Y, Gd, Pr, 
Yb) can be verified by experiment in the future. 

4.3 Specific heat and Entropy 

Once the phonon spectrum over the entire Brillouin zone is available, the vibrational heat capacity at constant volume 
Cy h can be calculated by using Eq. ([5]), while the electronic contribution to heat capacity at constant volume Cy can be 
obtained from the electronic DOS by using Eq. ([6]). Then, the specific heat at constant pressure C p can be computed by 
Eq. ([9]). As a comparison, both C p values, including electronic contribution and not, are plotted. We display these results 
in Figures [6] At high temperature, C y lh tends to the classical constant value 6R (R is molar gas constant), and Cy and C p 
still increase. Considering thermal electronic contributions to specific heat, we find that Cy lets C p suffer a little shift and 
can be negligible. This character can be understood from the electronic DOS at the Fermi level n(sf). AIRE intermetallic s 
have low electronic DOS near the Fermi level and the Fermi levels occur at a valley in the curves of electronic DOS |6|. 
Obviously, Cy is much smaller than the value C p - Cy = a 2 BVT. Hence, the electronic excitations affecting the thermal 
properties is inessential. Comparing with our previous study of MgRE intermetallic s l28ll , the electronic specific heat can 



7 



not be negligible since the Fermi level for MgRE intermetallics occur above a peak in the electronic DOS B31 and the 
electronic excitations affecting the thermal properties is remarkable. Whether we consider the electronic contribution or 
not, the theoretical C p keeps positive slopes obviously. At low temperature below ~ 200^, the discrepancy between C p 
and Cy can be neglected, and Cy values increase rapidly. We show temperature dependence of Cy for A1RE(RE= Y, Gd, Pr, 
and Yb) below 200K as in Figure [7] At low temperature, the specific heats Cy of four AIRE intermetallics obey the law of 
T 3 . The values of Cy are found A1Y< AlGd< AlPr< AlYb, and this character can be understood from elementary excitation 
of phonon relating to the mass of unit cell. We have the order of the mass of unit cell Maiy < ^AiGd < ^AiPr < ^AiYb- 
At low temperature limitation T — > OK, the specific heat Cy is determined from the vibrations of the low frequencies. 

The entropy as a function of temperature is an important thermodynamic quantity in thermodynamic modeling. The 
calculation results of entropies for the fore B2-type structures AIRE are shown in Figure [8] It can be clearly seen that 
the entropies of the four compounds increase with temperature. Through the temperature range - 1200^, the overall 
observation is that the order of the entropy S AlY<AlGd<AlPr<AlYb. At the temperature of 300^, the calculated 
values of entropy for A1Y, AlGd, AlPr, and AlYb are 67.48 JmoHK" 1 , 74.87 JmoHK" 1 , 77.93 JmoHK" 1 , and 85.03 
JmoHK" 1 . 

5 Conclusions 

In conclusion, the phonon and thermodynamic properties, such as the thermal expansions, the heat capacities at 
constant volume and constant pressure, the isothermal bulk modulus, and the entropy as function of temperature, of 
rare-earth-aluminum intermetallics AIRE (RE=Y, Gd, Pr, Yb) with B2-type structure are investigated by DFT and DFPT 
within the QHA. The phonon spectra and phonon density of states, including the phonon partial density of states and 
total density of states, have been discussed. Our results demonstrate that the density of states are mostly composed of Al 
states at the high frequency and composed RE (RE=Y, Gd, Pr, and Yb) states at low frequency. Through the calculated 
temperature range - 1200^, A1Y is the most incompressible and AlYb is the most compressible. The differences of the 
isothermal bulk moduli among the four intermetallics almost remain unchanged with increasing temperature. The thermal 
expansions of the four compounds are found AlYb>AlY>AlGd>AlPr, and those of A1Y and AlGd are nearly equivalent. 
The specific heats and entropy increase with temperature, while C y lb tends to a constant and C p keeps positive slopes at 
hight temperature. The electronic contribution to the specific heat is discussed, and the presented results show that the 
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thermal electronic excitation affecting the thermal properties is inessential. 
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Table 1: The equilibrium lattice constants ao, the isohtermal bulk modulus Bo, and the pressure derivatives of 
the isothermal bulk modulus B f at T = OK for AIRE (RE=Y, Gd, Pr, Yb) in our calculation in comparison with 
the previous calculated results and the experiment. 
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Figure 1: (Color online) The Helmholtz free energy F(V,T) as a function of unit-cell volume of A1Y. The 
circles denote F(V, T) at the volume points at every 100K between and 1200K. The solid curves show the 
fitted thermodynamic functions from Vinet EOS. The minimum values of the fitted thermodynamic functions 
at temperatures are depicted by the crosses. The dashed curve passing through the crosses is guide to the eye. 
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Figure 2: (Color online) Phonon-dispersion curves of (a) A1Y, (b) AlGd, (c) AlPr, and (d) AlYb. 
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Figure 3: Phonon TDOS and PDOS of (a) A1Y, (b) AlGd, (c) AlPr, and (d) AlYb. The solid curves, dashed 
curves, and dotted curves represent the phonon total DOS (TDOS), partial DOS (PDOS) of RE (RE=Y, Gd, Pr, 
and Yb) states, and PDOS of Al states, respectively. The PDOS indicates that the DOS are mostly composed 
of Al states at high frequency and RE states at low frequency. 
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Figure 4: (Color online) Isothermal bulk moduli B as a function of temperature. 
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Figure 5: (Color online) (a) Temperature dependence of the linear thermal expansion e(T) for A1RE(RE=Y, 
Gd, Pr, and Yb); (b) The coefficients of volume thermal expansion a for A1RE(RE=Y, Gd, Pr, and Yb) as a 
function of temperature. 
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Figure 6: Temperature dependence of heat capacity of (a) A1Y, (b) AlGd, (c) AlPr, and (d) AlYb. Solid and 
dashed curves denote the calculated C p , including electronic contribution and not, respectively. Doted and 
dot-dashed curves show vibrational and electronic Cy, respectively. 
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